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ABSTRACT 

A methodology for estimating physical parameters in a class of structural acoustic sys- 
tems is presented. The general model under consideration consists of an interior cavity which 
is separated from an exterior noise source by an enclosing elastic structure. Piezoceramic 
patches are bonded to or embedded in the structure; these can be used both as actuators and 
sensors in applications ranging from the control of interior noise levels to the determination of 
structural flaws through nondestructive evaluation techniques. The presence and excitation 
of the patches, however, changes the geometry and material properties of the structure as 
well as involves unknown patch parameters, thus necessitating the development of parameter 
estimation techniques which are applicable in this coupled setting. In developing a framework 
for approximation, parameter estimation and implementation, strong consideration is given 
to the fact that the input operator is unbonded due to the discrete nature of the patches. 
Moreover, the model is weakly nonlinear as a result of the coupling mechanism between the 
structural vibrations and the interior acoustic dynamics. Within this context, an illustrating 
model is given, well-posedness and approximation results are discussed and an applicable pa- 
rameter estimation methodology is presented. The scheme is then illustrated through several 
numerical examples with simulations modeling a variety of commonly used structural acoustic 
techniques for system excitation and data collection. 
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grants A FOSR-9O-0O9 1 and AFOSR-F49620-93-1-0198. This research was also supported by the National 
Aeronautics and Space Administration under NASA Contract Number NAS1-19480 while H.T.B. and R.C.S. 
were visiting scientists at the Institute for Computer Applications in Science and Engineering (ICASE), 
NASA Langley Research (-enter, Hampton, VA 23081. 
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1 Introduction 


The recent success of piezoceramic materials as sensors and actuators in applications involv- 
ing structural vibrations lias spawned intense study into questions regarding the modeling 
of piezoceramic actuator/sensor interactions with underlying structures, the optimal design, 
placement and number of actuators/sensors to be used, and the development of effective con- 
trol strategies in a variety of environments. When bonded to or embedded in a thin structure 
(beam, plate or shell), piezoceramic patches derive their actuating capacity from the piez- 
electric property that an induced voltage produces a strain in the material thus leading to 
the potential for producing in-plane forces and/or moments when the patches are mounted in 
pairs. Conversely, their sensing capabilities are due to the dual piezoelectric effect; namely, a 
mechanical force leads to the generation of a proportional voltage across the element which 
can then be used to measure accumulated strain. 

However, the bonding or embedding of patches in the underlying structure changes not 
only the geometry of the structure but also physical properties such as the density, stiffness, 
Poisson ratio and damping coefficients, when indeed, many of these parameters are unknown 
even for the homogeneous, uniform host structure material. This necessitates the development 
of effective parameter identification methods to he used when estimating system parameters in 
applications such as experimental model validation, the determination of optimal placement, 
and number of patches, the use of piezoceramic patches in nondestructive evaluation (NI)ft) 
techniques, and the implementation of model-based control schemes. We point out that the 
estimation of physical parameters in this setting differs from that considered in much of the 
previous literature (see [6] and the references therein) in that here the patch contributions 
to the system lead to unbounded (discontinuous) input and output operators due to the fact 
that the patches cover only discrete portions of the structure. 

Parameter estimation methods for distributed parameter systems involving unbounded 
operators have been developed and tested in the case in which piezoceramic patches are used 
as sensors and actuators when bonded to a transversely vibrating beam [15, 16]. There, fit- 
to-data techniques involving PI.) ft models were developed which could be used to estimate 
unknown beam parameters given various data, forms. Moreover, in that setting, results per- 
taining to convergence and continuous dependence oil data were obtained in a variety of cases 
involving physically tractable methods for exciting the system and measuring data. 

In this work, we develop an analogous methodology which can be used for estimating 
physical parameters in structural acoustic systems. In the systems of interest, an exterior 
noise source is separated from an interior cavity by a thin elastic structure (a beam, plate 
or shell). As energy is transferred from the exterior field to the structure, vibrations develop 
which then lead to unwanted interior noise through acoustic/structure interactions. Control of 
this unwanted interior noise is accomplished through sensing and actuating via, piezoceramic 
patches which are bonded to the structure. Before model-based control schemes can be imple- 
mented, however, the physical parameters of the structure (which now includes the patches) 
must be estimated from data which is collected both on the structure and from the acoustic 
response in the enclosed cavity. Although similarities exist between the problem of estimat- 
ing physical parameters for the isolated structure and that involving the structural acoustic 
system, the hyperbolic contributions due to the acoustic component sufficiently change the 



problem dynamics so as to warrant in-depth study of techniques for the latter coupled system. 
Finally, we note that although the initial impetus for developing distributed parameter esti- 
mation techniques for structural acoustic systems was motivated by model-based noise control 
considerations, the same techniques can be used when performing vibration analysis or using 
NDE methods to determine structural flaws in these coupled systems. 

The structural acoustic problem used here to motivate and illustrate the development of 
an appropriate parameter estimation methodology consists of a 2-D enclosed cavity which is 
separated from the perturbing exterior noise field by a thin beam. This model represents 
a 2-D slice from a 3-D model for several experimental apparata being used in the Acoustics 
Division, NASA Langley Research Center, to test modeling, parameter estimation and control 
strategies. We add, however, that the methodology being presented is equally valid for es- 
timating parameters in many 3-D models representing various experimental setups currently 
in use (see the models in [9, 10]). This 2-D model was chosen simply because it simplifies 
the discussion and more clearly illustrates the process involved in developing the parameter 
identification techniques. 

The model being used to illustrate the methodology is weakly nonlinear due to the manner 
in which the structural vibrations couple with the interior acoustic fields. While linearization 
provides a very good approximation to the system dynamics (see [1]), we retain the non- 
linearity here so as to illustrate some of the general analytic assumptions which are made 
when extending well-posedness and parameter convergence results for the corresponding lin- 
ear problem to a weakly nonlinear problem of this type. This also facilitates the demonstration 
of numerical techniques which can be used when simulating, testing and implementing the 
parameter estimation scheme in the nonlinear problem. In discussing parameter estimation 
methods for structural acoustics problems of this type, our emphasis is on the formulation of 
the problem in a manner which is conducive to approximation and implementation both in 
the linear and nonlinear forms as well as under a variety of damping assumptions. 

In the second section of this presentation, a model for the 2-D system being used to 
illustrate the parameter estimation method is presented. Details regarding the modeling of 
the acoustic and structural components as well as coupling conditions are given, and care is 
taken to motivate the assumptions which lead to various damping conditions in the system 
model. In formulating the strong form of the system model, details are also given regarding 
the interactions between the piezoceramic patches and the underlying structure (beam) as 
well as the weakly nonlinear interactions between the beam and the interior acoustic field. To 
provide a formulation which is conducive to approximation in the context of unbounded input 
operators as well as facilitates the discussion of well-posedness results, the weak form of the 
system equations is then developed and posed in terms of sesquilinear forms and the bounded 
operators which they define. Finally, we show that the solution trajectories can be expressed 
in terms of a semigroup on an appropriate space and within this framework, the assumptions 
underlying the well-posedness results for the linear and nonlinear problems are discussed. 

A discretization method suitable for simulations and the implementation of the parameter 
estimation method is outlined in Section 3. This discussion is kept relatively brief since details 
regarding the corresponding finite dimensional system for the linearized problem can be found 
in [1]. However, the section does provide additional details concerning the discretization of the 
nonlinear component of the operator and a brief algorithm for carrying out this discretizetion 
is included. 
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A parameter estimation scheme suitable for data consisting of displacement, velocity or 
acceleration measurements on the beam, voltage measurements from the patches, or pressure 
measurement inside the cavity is presented in Section 4. Assumptions on the form of the 
unknown parameters are discussed and conditions leading to convergence and continuous 
dependence on data results for the linear and nonlinear problems are outlined. 

In Section 5, numerical examples demonstrating the parameter estimation techniques with 
a variety of data types and methods for exciting the system are presented. Specifically, 
examples are given in which the force to the system is provided by a numerically simulated 
acoustic source, simulated voltage inputs to the patches, and a simulated hammer impact 
to the beam. The simulated data under comparison in the examples consists of acceleration 
values of the beam, voltage (accumulated strain) values from the patches, and voltage values 
in conjunction with interior acoustic pressure values. The conclusions from this study and 
physical considerations concerning the implementation of the method are summarized in the 
final section of the paper. 


2 System Model 

The model of interest consists of an exterior noise source which is separated from an interior 
cavity 0(f) by a common elastic boundary Po(f) that is modeled by an Euler-Bernoulli beam 
as depicted in Figure 1. The beam is assumed to have length (i, width b and thickness h. 
The Young’s modulus, mass density (in mass per unit volume) and damping coefficient for 
the homogeneous beam are denoted by Eb,pb and c Dbi respectively. Due to the nature of the 
exterior forces and the manner in which the patches are excited, we will be considering only- 
transverse vibrations iv(t,x). 

Bonded to the beam are piezoceramic patches which are mounted in pairs as depicted 
in Figure 1. In this discussion, it is assumed that the patches have thickness T , width 6, 
Young’s modulus E pe , density p pe , and damping coefficient c. Dpe . Moreover, it is assumed 
that the bonding layers for each patch have equal thicknesses, Young’s moduli, densities and 
damping coefficients, and these parameters are denoted by 7’^, Eb?, py and c-Dy, respectively. 
We emphasize that these assumptions are made solely for clarity of presentation, and similar 
results can be obtained in an analogous manner for the more general case in which the patches 
and bonding layers have differing thicknesses and material properties (see, for example, [12]). 

Finally, we assume that inside the cavity, there is a region 0 — Ul=u which provides 
a first approximation to the interior objects (eg., passengers, seats, storage compartments) 
which disrupt and damp the interior acoustic fields. This region is assumed to have positive 
measure and is taken to be small as compared to 0(f) (see Figure 1 where Nr — 4). 

2.1 Acoustic Component 

We consider first the acoustic wave dynamics inside the cavity 0(f). The variables of interest 
to us are the pressure P, density p, and the velocity U, each of which can be represented in 




Fig. 1. The 2- D domain. 

terms of a mean and fluctuating component 

P{t,x) = P 0 (x) + p( i . x) 
p(t,x) = po{x) + p(t,x) 

(/(t,,x) = U 0 (x) + ?<(/,.?) 

(here we have taken x = (x,j/) and are assuming that the rate of sound travel is sufficiently 
rapid so that little heat transfer takes place). For the range of magnitudes involved in these 
problems (< f 40 dB), it is usual to assume linear relations when considering constitutive laws 
and force balancing [20], and we make that assumption throughout the analysis which follows. 

In the region fi(2)/0, air damping is omitted due to the relatively small dimensions of the 
type of experimental cavities being simulated. Hence in that region, an increase in pressure 
brings about a proportional strain with the ratio defined as the bulk modulus of elasticity Ej 
where 

E f = -~ or p=-EjV-s (2.1) 

V • S 

(s(t,x) denotes the displacement of the center of gravity of an infinitesimal element of the 
medium and satisfies s t — u) . On the other hand, the material objects lumped in ft will 
provide some medium damping and here we assume that a change in pressure yields 

p — —EV - s — dS7 • s t (2.2) 

where E and el denote the bulk modulus of elasticity and damping coefficient of the medium 
in ft . We point out that this use of a generalized Hooke’s law in which stress is proportional to 



a linear combination of strain and strain rate is done under the assumption that relatively low 
acoustic frequencies are excited (< 1000 Hz), and is similar to the constitutive law leading to 
Kelvin- Voigt damping of vibrations in elastic, materials. We also emphasize that this damping 
model should be considered as a first approximation to the actual acoustic damping mechanism 
in the medium S 2 , and depending on the specific materials involved, the manner of acoustic 
excitation, and the geometry of the physical system, more comprehensive models may be 
required to accurately describe the medium damping. 

Force balancing in the acoustic cavity yields the relation (the Luler equation) 
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( 2 . 3 ) 
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denotes the equilibrium density of the medium ((2.3) is equivalent to the linearized momentum 
equation in fluid dynamics). By taking two time derivatives of (2.1) and (2.2), the divergence 
of (2.3), and eliminating cross terms, one arrives at 
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Taking tin* curl of the momentum equation (2.3) and noting that u = s t we obtain ^(V x 
p Q u) = 0. Hence the vorticity cD = V x p 0 u is constant in time in f 1(f). Under the assumption 
that tin* initial vorticity u>(0) is zero, we may conclude that V x p 0 u = 0 for all time or that 
the flow is irrotational in f2(/,). I bus in fi(/), there exists a scalar velocity potential O such 
that 

poll = —S7<p . (2.5) 


The flow is more complex near and in the region as a result of the viscous effects and 
medium damping. This can potentially lead to rotational components in the acoustic field 
which m general necessitates the use of a, vector potential. As a first approximation, however, 
we are assuming that the rotational components near and in Q are negligible and a relation 
of the form (2.5) will be used throughout the acoustic cavity U(/). 

Substituting (2.5) into the momentum equation (2.3) we find 


V 



= 0, 


which implies that ]> — since no acoustic sources or sinks are present in $2(/). Use of this 
pressure expression (actually p t = after differentiation in time once in equations (2.1). 
(2.2) then yields 


. — _\o 
W PS 


xe tt(t)/n 



Finally, hardwall boundary conditions are assumed since the experimental apparatus being 
modeled has concrete walls. With n denoting the outward unit normal to the cavity and c 2 (x) 
and d(x) denoting the speed of sound and damping coefficients given by 


c 2 (x) = 


E f/Pf 

E/p 


x e 

x £ 


the acoustic response can be modeled as 


d 2 <t> 

dt 2 


— c 2 A(j> “I - dA(f> t 


d{x) = 


0 

d/p 


x € 

x g , 


(x,y) £ n(t.) ,t> 0 , 


V (j) • n = 0 


(x,y) e r > o . 


( 2 . 6 ) 


We emphasize that this model was derived under the assumption that the only acoustic 
damping in the cavity occurs in the region $7 and hence d(x) = 0 in the rest of the chamber, 
i.e. for x £ fi(7)/17. Moreover, we have assumed that the flow is irrotational in the region 

n(f). 


2.2 Beam Component 

Through force and moment balancing, the equation of motion for the transverse displacements 
w of the beam are found to be 


d 2 M 


pw,,+ dx 2 

ru(f,0) = 0) = w(t,a ) 


(t,x) ~ -<f> t (t,x,w(t,x)) + f{t,x) 

div 

dx v-,~> fa 


xO < x < a , t > 0 , 
(t,a) — 0 t > 0 , 


where M is the total beam moment, / is the force due to pressure from the exterior noise field 
and x, w(t , x)) is the backpressure due to the ensuing acoustic waves inside the cavity 
(this latter term is in general nonlinear since its effect occurs on the surface of the vibrating 
beam). For pairs of patches having edges x i and x 2 , the density of the structure is 


p(x) pbhb “f~ 26 ( pbfTbf ~b Ppe'E') yp e (.c) 
where the characteristic function 


(2.7) 


Xpe(x) = 


1 

0 


X, < X < X-2 
otherwise 


localizes the patch effects between the endpoints X] and x 2 (see [11, 12] for details). 
The general beam moment 


M. (t, x) = M ( t , x) — Mpe(t , x) 


consists of an internal component M, depending on material and geometric properties of the 
beam and patches, and an external component M pe (the control term) which results from the 
activation of the patches through an applied voltage (see Figure 2). 
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For a beam undergoing pure bending motion with out-of-phase excitation of the patches, 
the internal and external moments are given by 




EI(x) 


8 2 w 

~dx> 


+ C D I(x) 


d :i w 

dx 2 dt 


M pe (t,x) = k' B VXpr(x) 


where V is the voltage into the patches. As shown in [1 1], the stiffness, damping and control 
constants for the combined structure are 

E I (xj Efo y ~ "F ~F Ep e (l 3pe ] Xpe(*r) 

, h^b 2 b r i 

c D l{x) - + — [CDbt.asbf + c D p e a 3pe j x P e(*) 


(2.8) 


K b = E p M:n(h + 27 ’ m + T) . 


The constants a^y and a^ pe given by — (h /2 + Ty) 3 — (/?./ 2) 3 and a 3pe — (h/ 2 + Ty + T) 3 — 
(h /2 + Ty ) 3 result from the integration of stresses through the bonding layer and patch, and 
f/31 is a piezoceramic constant which relates the amount of strain produced in the patch to 
the level of voltage being applied. 

At this point it is worth commenting further on the damping term co pe which is taken to 
be a combination of the Kelvin-Voigt. damping coefficient for the patch and the damping which 
results from the production of current when the structure vibrates. This latter contribution to 
the damping results from the piezoelectric effect of the patches which dictates that a voltage 
is produced when the patch is subjected to in-plane strains. Under the assumption that the 
Kelvin-Voigt (material) and electrical damping have approximately the same types of effect 
in the patch, we have combined the two into the coefficient cp pe which must be considered to 
be unknown and like the other parameters, must ultimately be estimated using data fitting 
techniques with experimental data when considering actual applications. We also point out 
that the expression (2.8) can easily be generalized to include the possibility of differing material 
properties in the two patches or bonding layers (again, see [11, 12]). 

The fact that the patches generate a voltage when strained implies that they can be used 
as sensors as well as actuators. Specifically, the voltage produced by the patches when the 
beam bends is 


V s (t) = K. 


Jxi 


n 2 d 1 


W 


dx 1 


-(/, x)dx = K: 


dw . , dtv , 


( 2 . 9 ) 


where the sensor constant /C 5 depends upon piezoelectric material properties as well as the 
geometry and size of the patch (see [ 18 ]). Hence the voltage provides a measure of accumulated 
strain in the beam (see Figure 2 ) thus enabling the patch to serve as a sensor in a variety of 
applications involving the measurement of beam vibrations. 
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Fig. 2. Acoustic cavity with piezoceramic patches; (a) fixed domain it and boundary I |J To; (b) patch 
generation of pure bending moments , and (c) patch sensing of accumulated strain . 


2.3 Coupling Conditions 

In the model discussed thus far, the structural and internal acoustic responses are coupled 
through the backpressure <f>t(t, x, x)) on the surface of the beam. A second coupling 
equation is the continuity of velocity (or momentum) condition 

w t (t , x) = — V </>(/, x, w(t, x)) • j , 0 < x < a , t > 0 

Pf 

which results from the assumption that the beam is impenetrable to air. We point out that the 
velocity condition provides a form of damping to the beam which is similar to that obtained 
with the incorporation of viscous (air) damping effects (modeled by a term of the form 7 w t in 
the beam equation). As noted in the examples, the internal Kelvin- Voigt damping in concert 
with the coupling effects due to the continuity of velocity and backpressure causes a beam 
response which differs somewhat from that observed with an uncoupled, undamped beam 
having the same dimensions. 

The model we have developed has nonlinearities in the (i) variable domain (ii) back 

pressure term — x, w(t, x)), and (iii) velocity term ~ V^(/, x, iu(i, x)) ■ j. Under an as- 
sumption of small displacements (w(t,x) = w(t,x) + S where ih = 0 ) which is inherent in 
the Euler- Bernoulli formulation^ the variable domain f 2 (/) can be approximated by the fixed 
domain Q = [ 0 ,a] x [ 0 , f] as shown in Figure 2. Note that with this assumption, the velocity 
term — x, w(t, x)) * j can be approximated by the normal term V<p(t,x , w{t.,x)) ■ 11 

which arises when developing the weak form of the equation. The fully nonlinear form of the 
back pressure coupling term is retained throughout the following discussion. 
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2.4 Strong Form of the System Equations 

For the coupled system in which s pairs of patches are bonded to the beam and excited out- 
of-phase, the acoustic, structural and coupling components just discussed can be combined to 
yield the approximate nonlinear model 


<j> tt — c 2 A<j> + dA(f> t (;r, y) G J7 , / > 0 , 
V <f> ■ h = 0 ( x ,;// )€ F , / > 0 , 


— V(p(t, , x, u)(t, x)) • h = w t (t , x ) 0 < x < a , t > 0 , 

Pi 


d 2 

+ o? 


El 


, () 2 w 


dx 2 


+ ('i)l 


(Pw \ 

dx 2 (H ) 


= -od/.-r- w(t ,./•)) + /(/.,; r) + ~ u t (t)x pf ,{x) 

dx ,=i 

w(t.O) = 0) = ^l!(t,a) = ——(/,«) = 0 t > 0 , 

ax Ox 

<t>{0,x,y) = <p 0 {x,y) , w(Q,x) = w 0 (x) 

(pt{0,x,y) = <j) : {x,y) , iO((0,;r) = ?e,(.r) . 


( 2 . 10 ) 

0 < x < a , 
t > 0 , 


Here u t (t) is the voltage being applied to the i th patch and \ v( t denotes the characteristic 
function over the i ih patch. 

We point out that the piezoceramic material parameters K.f (see (2.9)) and K f , i = 
as well as the beam parameters p,El and cpl arc considered to be unknown and 
are estimated using inverse problem techniques as discussed in later sections. While the 
expressions given in (2.7) and (2.8) can be used. as starting values in the parameter estimation 
routines, experimental evidence (see [15, 16]) has indicated that the final parameter values 
can vary quite significantly from the analytic values due to the contributions from the bonding 
layer, variation in the measurement of physical constants, and nonuniformities in the various 
materials. This combined with the lack of analytic expressions for the damping constant 
necessitates the estimation of these parameters before model-based control strategies can be 
implemented. 

We also emphasize that the parameters Kf,p,EI and cpl are piecewise constant in na- 
ture due to the presence and differing material properties of the bonding layer and patches 
(see (2.8) as well as the results in [15]). This leads to difficulties with the strong form of the 
system equations since it necessitates the second derivatives of the Heaviside function (equiv- 
alently, derivatives of the Dirac delta) thus yielding an unbounded control input operator. 
The differentiation of the discontinuous material parameters also leads to difficulties when ap- 
proximating the dynamics of the coupled system. To avoid these problems, it is advantageous 
to formulate the problem in weak or variational form (the use of the variational form also 
permits the use of basis functions having less smoothness than required for those used when 
approximating the solution to the strong form of the equations). 
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Finally, we note that in the case of no acoustic cavity damping (d = 0), the model (2.10) 
is completely equivalent to the nonlinear models that are the basis of the investigations in 
[1, 8, 10, 11] if one replaces the coupling terms pj<j) t and v = — V0 in those models by (j>t and 
v = — respectively. That is, the potential used in those references differs from the one 

used here by a multiplicative factor 


2.5 Weak Form of the System Equations 

An appropriate choice for the state of the second-order problem (2.10) is the pair (0, u>) 
consisting of the acoustic potential and beam displacement. It follows immediately that with 
this choice, the state for the problem in first-order form is Z — (</>, w, 0, te) which contains the 
pressure (since p = <j>) as well as the beam displacement and velocity. 

The state space and space of test functions are taken to be the product spaces 7i — V x H 
and V = V X V where the Hilbert spaces H and V are given by H = L 2 (0) x L 2 (Fo) and V = 
H l {St) x //q(Fo). Here L 2 (il) and H*( fl) denote the quotient spaces of L 2 and Z/ 1 over the con- 
stant functions and //^(To) is given by Hq(F 0 ) = {0 £ H 2 (To) : — VA^) = 0 at x — 0,«} 

(the use of the quotient space results from the fact that the potentials are determined only 
up to a constant). From energy considerations, the // and V inner products are taken to be 

((») ’(»)) „ = + 

/ f El D 2 wD 2 r)d i y 

Jn J r 0 

while the product space inner products are given by 

((*)'( A )) h = <4M0, + <*,A)„ 



Integration by parts in a manner analogous to that in [1] in combination with the approx- 
imation of the variable domain fi(£) by the fixed domain il = [0,a] X [0J] and the use of 
Green’s theorem then yields the nonlinear first-order variational form 


+ / V<f) ■ + f El D 2 wD 2 i)d~f 

Ju J r 0 

+J -d\ ? (/)-V('du>-\- j {c d ID 2 wD 2 t) + </j ( a;)?/ — ti^J | d~f 
= f i2 K f u '( t )Xpe l D 2 7]dj + f fl)d-j 


(2.i i ; 
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for all ^ = (£,r/) G V 7 . We point out that in this variational form the derivatives have been 
transferred from the plate and patch moments onto the test functions. This eliminates the 
problem of having to approximate the derivatives of the characteristic function and the Dirac 
delta as is the case with the strong form of the equations. 

The system (2.11) can be formally approximated by replacing the state variables by their 
finite dimensional approximations and constructing the resulting matrix system, and it is in 
this form that we will consider approximation strategies in Section 4. In order to discuss the 
well-posedness of the problem, however, it is advantageous to pose the problem in an abstract 
Cauchy format as discussed in the next subsection. 


2,6 Abstract First Order Formulation 

Following the theoretical work in [2, 3], it is advantageous to formulate the problem in terms 
of sesqui linear forms and the bounded operators which they define (see also [1, 10] for further 
examples detailing the abstract formulation of structural acoustic systems in this manner). 
We begin by pointing out that the Hilbert spaces H and V form a (Jelfand triple V c * H ~ 
H* <— > V * with pivot space // (further details concerning the basic definitions and fundamental 
functional analysis theory here can be found in [22]). For <J> - and = (£,//) in V\ we 

then define the sesquilinear forms cr t : V X V — ► C, i = 1,2 by 


o - ! (4>, = / V (f) • V £ (Ilo 4” [ EID' 2 wD 2 7](h , 

J 0 *7 To 

< 7 - 2 ( 0 , '!') = / — 2 <i V(f) ■ + f {c n ID 2 wD 2 Tf + \<f>T} - »’£]} d-y 

J 12 C J Fq 


(2.12) 


With these definitions, it is straightforward to show that (7\ and cr-i satisfy various conti- 
nuity, symmetry and coercivity conditions. Namely, satisfies 


|(7](4>,^)| < c t |<I>|v'|'P|v , for some C] G IR 
Re cti( 4>,$) > c-il^lv , for some c 2 > 0 


= cr, (vp, <I>) 


(bounded) 

( V / -elliptic) 
(symmetric) 


(2.13) 


for all <K vj/ g V (the boundedness results from Schwarz’s inequality for inner products in 
conjunction with equivalence results for various Sobolev norms, while the K-ellipticity and 
symmetry of ctj follow directly from the fact that (r,(4>, 'P) = (4>, V P) V ). 

The behaviour of <r 2 depends upon the form of damping in the cavity fl. If damping is 
included in a region fl having positive measure as discussed in Section 2.1, then <r 2 satisfies 


|<7 2 (4>, T)| < C;j | <1* | v/ 1 T | v , for some c 3 G IR (bounded) 

Re cr 2 (4>, 4>) > c *4 14> |^r , for some c 4 > 0 ( V-elliptic) . 


(2.14) 
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The C-ellipticity follows from the fact that 


Re <r 2 ($, 4>) = [ I Vf\ 2 du + 

Jn c l 

> k x [ \V<j>\ 2 du + k 2 f El (D 2 wYd 1 

> C A 

with the final inequality resulting from the observation that there exists a constant ky, such 
that k i /jj |V<^| 2 (fo; > Jq \\7<j)\ 2 du) for all 4> € V. On the other hand, if the cavity is taken 
to be completely open with no medium damping assumed (hence meu.s(fi) = 0 and d = 0 
throughout f 2), then one can only establish the weaker //-semiellipticity condition 

Re > c 4 \$\' 2 H , c 4 > 0 

for (72 (in this case c 4 = 0 since contains no acoustic components when damping is 

omitted in the cavity). As discussed in the next section, the inclusion of medium damping 
within the subregions (1 and the resulting K-elliptic behavior of fj <2 yields stronger semigroup 
results than those obtained when acoustic damping is omitted. 

To account for the patch contributions when a voltage is applied, we let V denote the 
Hilbert space containing the voltage inputs and we define the operator B £ £((/, V *) by 

{Bu, H>) v . v = f K.f UiXpei 

Jr ° <= i 

for 'P £ V 7 , where is the usual duality pairing. Finally, the external forcing and 

nonlinear perturbation terms are given by F = {0,f/pb) and G(z,z t ) = (0,— </>t{w)) where 
<f>t{w) = (j>t{t,x,w(t,x)) = <f> t (t,x,w(t,x)) — (j)t(t,; r,0) denotes the nonlinear perturbation to 
the linear coupling term. 

With these definitions, we can write the system (2.11) in the abstract weak or variational 
form 


j r c D l(D 1 wj 2 d 1 


(zu(t), *) V .' V + <rt(z t {t), #) + <T, (*(/), ' 0 = (Bu(t) + F + G(z , z t ), V) v . y (2.15) 

for in V. We reiterate that the state for the problem in second-order form is given by 
z(t.) - (<f>(t, •, •), w(t, •)) in V «->■ H. 

To pose the system in first-order form, we form the product space terms Bu(t) — (0, £?«(/)), 
T{t) = (0 ,F(t)) and G{Z(t)) = (0, <7(ar(i), ^t(<))) in V* = C x V* and define the operators 
c £(V,V*) by 

for i — 1,2 (the existence of A\ and A 2 is guaranteed by the boundedness of and a 2 ). 
Then, for the state Z(t) = (z(t. ), 2 ( (/)) = {(j) , w , 6 , in) in 7i = V x //, the weak form (2.15) is 
formally equivalent to the system 

Z t (t) = AZ(t) + C(t.Z(t)) (2.16) 
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in V* where 


c(t,z(t)) = Bu(t) + ?(t) + g(Z(t)) 


(2.17) 


and 

dom A = {0 = (T, A) e H : A <E V, /1,T + ,d 2 .\ € H) 

0 / j (2-l«) 

/l] —A 2 

(further details concerning the formulation of the first-order system in the linear case can 
be found in [1]). The representation (2.16) is formal in the sense that the manner in which 
differentiation and the resulting solution exists has yet to be specified. This will be discussed 
next in the context of proving well-posedness results for the system model. 

Finally, in discussing well-posedness results and parameter convergence in the sections 
which follow, it proves useful to compare the nonlinear system with that which results when 
both coupling terms are linearized. The latter is found by replacing the term <j) t (t,x, w) by 
<f)t(t,x,0) in the various expressions for the coupled system (equivalently, take (i(z,z t ) = 0). 
This yields the second-order variational form 

MO, # W + *) + *i(*(0, #) = + /\ *) v . y (2.19) 

for ^ £ V and consequently the first-order system 

Z t (t) = AZ(t) + Bu(t) + F(l) (2.20) 

in V*. 

2.7 Model Well-Posedness 

The first step in determining the well-posedness of the semilinear system model is to argue 
that A generates a semigroup on 7 H. As noted earlier, the sesquilinear form <7\ is V-elliptic, 
continuous and symmetric while a 2 is continuous and V -elliptic if damping is included in the 
region Cl or //-semielliptic if damping is omitted in the cavity. In both cases, the Lumer- 
FMnlips theorem (with arguments similar to those given in pages 82-84 of [2]) can be used to 
prove that the operator A defined in (2.18) generates a Co-semigroup T on the state space 7 H. 
The semigroup satisfies the exponential bound |T(Z)| < t ,ujl for t > 0 (where in fact, lj = 0 due 
to the fact that A is dissipative as shown in [2]). Moreover, if medium damping is included 
in the region Cl (which implies that rr 2 is V-elliptic), the semigroup T is analytic on TL. 

For the problem thus posed, the state lies in 7i which implies that the semigroup T 
generated by A : dom A C 7Y — > 7Y is defined on 7 i. The nonlinear forcing term C(t,Z(t)) = 
Bu(t) + J r (t)-\-G(Z(t)), however, lies in V* rather than H since the control term B £ £(//, V *) 
defines the product space control term Bu(t) = (0, Bu(t)) £ {0} x V * C V x V* = V*. 
To remedy this, “extrapolation space” ideas and arguments similar to those' presented in 
[2, 4, 19] are used to extend the semigroup T(t) on TL to a semigroup T(t) on a larger s[>ace 
W* D {0} x V* so as to be compatible with the forcing term. 
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As detailed in [10], the space of interest is defined in terms of dom A* where 
dom A* = ix = ($,'!') e H\* e V, A]<i> - e H) 




-* ) 
A\*-A\* ) 


Specifically, the space VV = [dom *4*] is taken to be dom A* with the inner product 


= {(\ 0 - A*)*, (\ 0 -A*)*) n 


for some arbitrary but fixed A 0 with A 0 > cj (recall that the original solution semigroup 
satisfies the bound \T (t)\ < e wi ). As proven in [4], the resulting W norm is equivalent to 
the graph norm corresponding to A* . Moreover, we have that {0} X V* C W* = [dom .4*]* 
(see [10] for details). 

From the definition of A* and the equivalence of the W norm with the graph norm corre- 
sponding to A*, we can define *40 E VV* by 

(ie) W = (e,A)« 

for all 0 6 7-f, x £ W. With this definition and the Riesz representation theorem, it is 
shown in [10] that A is an extension of the original operator A from dom A C Ft to all of 77 
Finally, as proven in [4], the operator A is the infinitesimal generator of a Co-semigroup T(t) 
on VV* which is an extension of T(t) from Ft to VV* (note that T(t) is also analytic if medium 
damping is included within the region Cl). 

Having extended the operator A and hence the generated semigroup to a space which is 
compatible with the forcing function, we are now in a position to discuss criteria on the input 
terms F and B which guarantee the existence of a unique solution to the system model. In 
the corresponding linear problem, under reasonable regularity conditions on t >— » u(t) and 
t i — > F(t) } one can immediately argue the existence of a unique strong solution to the system 
in terms of the extended semigroup T ( t ). For the semilinear problem of interest, however, the 
nonlinear nonhomogeneous terms must satisfy certain continuity criteria in order to obtain 
similar results. For example, if we let X denote the reflexive Banach space W* and assume 
that C:[0,T]xX“)X defined in (2.17) is continuous in t on [0,7] and uniformly Lipschitz 
continuous on X , then the integral equation 


Z(t)=f{t)Z 0 + 

Jo 


0 

Bu(s) + F(s) + G(Z(s)) 


ds 


( 2 . 2 1 ; 


is well-defined for Bu + F + G{Z ) E L 2 ((0 ,7), V *). Moreover, for 2(0) = Z Q , the solution 
Z(t) of (2.21) is a unique mild solution to (2.16) (see Theorem 1.2, page 184 of [21]). In 
addition, if C : [0, 7] x X — ► X is Lipschitz continuous in both variables, then it follows from 
Theorem 1.6, page 189 of [21] that (2.21) provides the strong solution to (2.16) interpreted in 
the VV 3 " sense. 

The required continuity of the nonhomogeneous terms Bu and F is demonstrated in [10] 
and hence the remaining question concerns the Lipschitz continuity of the nonlinear coupling 
term G(z,z t ) — (0, — <j) t (w)). If we assume that the input terms F and Bu are sufficiently 
smooth so as to assure the necessary continuity in G(z,z t ), then the nonlinear system is 
well-posed. These results are summarized in the following result. 
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Lemma 1. (Well-Posedness of the Nonlinear System) Consider the nonlinear system 
represented by (2.15), (2.16) or (2.21) and assume that F and Bu are sufficiently smooth 
so that C(t,Z(t )) = Bu(t) + ^F(t) + Q(Z(t)) is Lipsrhitz continuous in both t and Z . 
Then (2.21) with Z^ £ 'H is the unique strong solution to (2.15) both when acoustic 
damping is assumed in the region fl and when no damping is present in the cavity (in 
which case, a 2 is only //-semielliptic). 

These results can be further extended when acoustic damping is assumed in the region fi 
and hence is V-elliptic and T is analytic. In this case, the mild (and thus strong) solutions 
to (2.20) are guaranteed to be equivalent to the weak or variational solutions to (2.19) for 
sufficiently smooth nonhomogeneous terms. The following theorem summarizes this result for 
the linearized problem. 

Theorem 1. (Linearized System: Equivalence of Solutions) Consider the system rep- 
resented by (2.19) or (2.20) and suppose that the mappings t ► u(t) and t >— > F(t) from 
[0 ,T] to IR 1 and V*, respectively, are Lipschitz continuous. Furthermore, assume that 
medium damping of the form (2.2) is present in the region Q C fl. Then for each 
Zo G Ft = dom A , we have that (2.20) taken with Z( 0) = Z 0 has a unique strong 
solution given by (2.21) with Q{Z) — 0. Moreover, this strong solution is equal to the 
weak solution of (2.19). 

The proof of the equivalency between strong and weak solutions follows that given in [3] 
for general second-order systems with unbounded input terms. We point out that numerous 
numerical results have indicated similar results for the nonlinear problem and the case in which 
damping is omitted in the cavity even though we have not extended the results equating the 
strong and weak solutions to cover those cases. 


3 State Approximation 


The modeling and well-posedness discussion thus far has been for the infinite dimensional non- 
linear structural acoustic system. In this section we discuss a Galerkin scheme for discretizing 
the problem which can be used when simulating the system dynamics, estimating the physical 
parameters, and calculating control gains (see [1]). This is accomplished by approximating 
the beam displacement and acoustic potential by spline and spectral expansions, respectively. 

3.1 System Approximation 

As detailed in [1] where the corresponding linear problem was considered, cubic splines are 
suitable for discretizing the beam displacement since they satisfy the smoothness requirement 
as well as being easily implemented when adapting to the fixed-end boundary conditions and 



patch discretizations. Specifically, the approximate beam displacement is taken to be the 
linear combination 

»*(«.*) = 2 "-ftOSTW 

1 = 1 

where B 71 is the i th cubic spline which has been modified to satisfy the boundary conditions. 
The acoustic potential is approximated by the Clalerkin expansion 


Wy 7/1 x 

j= 0 *=o 


where P?(x) and Pf(y ) denote the standard Legendre polynomials that have been scaled 
by transformation to the intervals [0, a] and [O,/ 7 ], respectively. The condition i + j ^ 0 
eliminates the constant function thus guaranteeing that the set of functions is suitable as a 
basis for the quotient space. We take { B k(x , O)}™^ — { P? (x)P-(y)} where m = 

(m x + l){m y + 1) - 1. 

The m and n — 1 dimensional approximating cavity and beam subspaces are taken to be 
H 7 c n = span {B'' 1 }™) and = span {£?”} "T, 1 , respectively, where B” 1 and Bf are the i th 
cavity and beam bases described above. Defining N - n — 1 , the approximating state 
space is then taken to be H N = //” l x H b and the product space for the first order system is 
U N = H N x H n . 

By restricting the infinite dimensional system (2.11) to 7i N x7i N , one obtains the nonlinear 
finite dimensional system 


with 


and 


' m, n 

0 ■ 


’ d N (t) ’ 


0 

M? . 


j) N (t) 


' A/f 

0 ' 


' t? N (0) ' 

r 

0 




L 




9? 

92 


-A? («."(!)) 


d N (t) 


0 

B n 


u(t) 


M? = diag[M", M"] , 

M? = diag[M A , M»\ , A ^ (>(*)) = 


— d i ag [ xt , Afy 


B N =[o,B?] T , F N {t) 




o ■ 




0 

F N [t) 


The vector ti N (t) = <f>%(t), • • • , w?{t),w?{t), • • • , {t)) T contains the N x 1 = 

(m + n — 1) x 1 approximate state coefficients while u(t) = • • • u s (t)) contains the s 

voltage values. The matrices M.£ and A^ are the mass and stiffness matrices which arise when 
solving the uncoupled wave equation with Neumann boundary conditions while M%f 2 ,Ay 2 and 
A 22 are the mass, stiffness and damping matrices which arise when solving the damped beam 
equation with fixed boundary conditions. The matrices M[{ and My 2 result from the choice 


16 


ot V inner product (see [1]). The contributions from the coupling terms are contained in the 
matrix and operator A^ 2 ( w ,,Si (/)) while the control, forcing and initial terms are contained 

in B 2 ' (t) and i)q , respectively. A more detail description of the various component 

matrices can be found in [1]. 

The equivalent finite dimensional Cauchy equation is then 

y N {t) = A N (y N {t)) + B N -u(t.) + F N (t) . | 

?/ v (0) = fjo 

with the various operator definitions following from those given above. As with the infinite 
dimensional problem, it will also be useful to consider the corresponding problem 

y N (t) = A N y N (t) + B N u(t) + F N (t) 

v v {:l2) 

2/ A (0) = .Vo 

when considering issues such as parameter convergence. The linear operator is obtained 
by considering the linearization of A^ 2 (w^ (t)) as discussed in [1], 


3.2 Algorithm for Constructing the Nonlinear Component 

The determination of state trajectories involves the repeated construction of the nonlinear 
operator (y^ (/ )). Although most of the components of this operator need to be constructed 

only once, the (n — 1) X m matrix A^ 2 ( w^ (f )) must be determined at each step in the' solution 
of the ordinary differential equation (ODE) (3.1) and hence its formation must bo made as 
efficient as possible. 

We first point out that for i,j = 1,* • • . m x ,m VT this (n — 1) x m matrix has the entries 


A 32 (« ;A 


hi - 1 


PAJ 


£ b;(x)P:'(x)p; wH(t)B'j(x)j dx 
where the integrals are evaluated via a Gaussian quadrature rule of the form 

[ f(x)dx « Yj c kf(*k) ■ 

J o rT 


Here n q is the number of quadrature points and c* is the quadrature weight. Those quadratures 
and hence the formation of the matrix can be efficiently accomplished in the manner outlined 
in the following algorithm. 


Numerical Algorithm for Creating the Nonlinear Component A£ 2 (/)): 

1. Create the matrices Ay la , A^ 2aw and A^ 2bj} which have the components 

- , [> 45 -.],, -<*«?(*<) . [A&hL - 


a n 

32a 


ik 


(hence dim A% 2a = n q x(n- 1 ) , dim ^^ 2au , = (n- 1 ) xn g and dim = n q x (?n J .+ 1 )). 

These matrices need to he formed only once and can be created before' solving the ODE 
system. 
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2. With the notation o denoting Hadamard or componentwise multiplication, perform the 


following operations in the ODE solver: 


* = A* 2 aW N (t) 
for j — 1 : niy + 1 

, dim x — n q X 1 

b = Pj{x) 

, dim b — n q X 1 

Am bj2 (w N (t)j = b * oncs( 1 , m x + 1) 

, dim Ay lhfl (ie w (/.)) = n q X (m x + 1 ) 

"4.326 j ( ti,Af (0) = ^32bj\ 0 ^32t >j2 { U,N (0) 

, dim Ay ib - (u) N (tj^j = n q x {m x + 1) 

ind 1 = (j - 1) • (m x + 1) + 1 
in d '2 = j ■ ( m r + 1 ) 

A 326 ( U,/V (0) (‘,indl : ind'2) = A$ 2bj (ui N (i 

!)) > 4im A% ib (t /(i)) = », X m 

end 


"4.32 (« jJV (0) = P! A 32aw ' A ^ 2 b (^C)) 

, dim A 32 ( u?jV (t)) = (n — 1) X m . 

By creating the matrices Af 2a , A$ 2aw and A 22h ^ offline, the time needed to solve the ODE 
system is reduced thus improving the efficiency of the scheme for parameter estimation and 
control applications. 


4 Parameter Estimation Problem 

The goal of the parameter estimation problem is to determine the “true” material parameters 
pb, EI,cpI and JCf , K'f pi = 1, • • • ,.s, given data measurements z from some observable sub- 
space Z of the state space. To pose this mathematically, we let q — ( pb , El , cpl , K.f , • • • , K.f , 
k'.s , • • • , K's ) and assume that q 6 Q where Q denotes an appropriately chosen admissible 
parameter space. The parameter estimation problem is to then seek q £ Q which minimizes 


•/(</) 


C 2 \Cx{Z{tuq)} 


(4.3) 


given pointwise temporal measurements 5, = z(ti) at given points on the beam and in the 
cavity. Note that this minimization is performed subject to Z = (<j),x,<f>,w) satisfying the 
coupled system equations (2.11) or (2.15). Depending on the experimental apparatus, the 
data observations may consist of position, velocity, acceleration, or accumulated strain mea- 
surements at points on the beam as well as pressure measurements inside the cavity. The form 
of the operator C- 2 depends upon whether one is performing the estimation procedures in the 
time domain or in the frequency domain. In the time domain, C 2 is the identity whereas it is 
the Fourier transform for estimation in the frequency domain. 

For time domain estimation with data consisting of position, velocity, or acceleration 
measurements at points x on the beam, the fit criterion functional to be minimized is 


-/(<?) = E 

l 


(Ew 

dt u 




(4.4) 


18 



with v — 0, 1 or 2. On the other hand, if a patch with endpoints at x x and x 2 is used to collect 
accumulated strain measurements, an appropriate functional is 


■H'l) = E 


K' s 


l 


( d w \ 


dw 

dx 




(4.5) 


(see (2.9)). Here the data consists of the voltage 1 measured across the patches. The patch data 
can also be combined with pressure measurements ]\ at points (.r. y ) in the cavity to provide 
a fit criterion functional 


^) = E - 


+ | </> t (ti,x,y) - pi 


4.6) 


which in some cases has more sensitivity than that in (4.5) which considers only strain mea- 
surements. 

The above fit, criteria can also be used with data that has been transformed to the frequency 
domain (in which case C 2 is the Fourier transform), and this is indeed a common procedure 
for starting the optimization process with data in which several frequencies are excited (see 
the comments in the examples as well as [13, 14]). In this case, optimization is qualitatively 
performed in the frequency domain until frequencies match sufficiently so that the optimization 
routines will converge with the time domain data. 

To facilitate the estimation of the material parameters p, El and coE we now make some 
assumptions regarding their spatial behavior. Because the beam and patches are considered 
to be homogeneous as well as uniform in width and thickness, it is reasonable to assume that 
the density, stiffness and damping parameters of the combined beam/piezoceramic patches are 
piecewise constant in nature (see for example, [15]). A suitable partition is then taken to be 
{xfc} = {0,a} U{ a i/ }}=i 2 5 where the 2,s points {o^} are the endpoints of the s piezoceramic 
patches. Finally, we assume that these parameters have the form 

p{x) 

EI{x) 
cdI(x) 


25+1 

= c kBk{x) , 'Cl = C 3 = • • • = c 2s + 1 

k = 1 

2.s+l 

= ^k&k{x) , Ci — C 3 = ‘ 1 * = <-2s+\ 

k= 1 

25+1 

= ^2 9 E = C 3 = • ‘ ' = C-2. s +l 

k=\ 


where the piecewise constant basis functions are defined l)y B k (x) = H(x — : r^_i) — H(x — 

The coefficient constraints c x = c 3 = * • • = c 2s + 1 , and so on, result from the uniformity of the 
beam in areas not covered by patches. Finally, we recall from the definitions (2.8) and (2.9) 
that the patch parameters /Cf , • * • , /Cf , and A f, • • * , A f are simply constants which depend 
on piezoelectric properties, the geometry and size of the patch, and patch and bonding layer 
properties. 

Although Q is finite dimensional with the above assumptions on the parameters, the min- 
imization of the fit criteria in (4.4), (4.5) and (4.6) involves an infinite dimensional state 
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and hence is not immediately tractable. With the Clalerkin schemes of the the last section, 
however, corresponding minimization problems involving the state approximations can be de- 
veloped and used when estimating physical parameters with these fit-to-data techniques. With 
w N , 4 and H N C V defined in the last section, finite dimensional functionals corresponding 
to those in (4.4), (4.5) and (4.6) are 


7) = E 

J N iv) = E 


d u w 


N 


at v 




, v = 0 , 1 , 2 , 


s (dw N dw N \ 


and 

= E 


y . s fdiv N , , dw N , \ „' 2 

* I — —(U,zuq) - Zi 


dx 


dx 


+ 


$ {U, x, y) - Pi 


(4.8) 

(4.9) 


(4.10) 


respectively. The approximate beam displacement w N and acoustic, pressure pf<j>? at the 
various points are found by solving either the nonlinear finite dimensional system ( 3 . 1 ) or the 
system (3.2) if one is considering the linearized problem. 

The following theorem taken from [16] specifies conditions under which convergence and 
continuous dependence (on data) of the solutions to the linearized finite dimensional parameter 
estimation problems involving the functionals (4.8), (4.9) and (4.10) can be expected. 


Theorem 2. (Linearized System) Let Q be a a compact subset of a metric space Q with 
metric d and assume that H N C V approximates V in the sense that for each <t> E V 7 , 
there exists E H N such that 


<J> 




V <^N) 


0 as N — > oo . 


(4.11) 


Furthermore, assume that cr^r/) and cr 2 ((/) defined in (2.12) are K-elliptic, continuous, 
and satisfy the continuity with respect to parameter condition 


|<Ti((/)(4>, 4») - (Ti(q)($, ty)| < 7 t d((/,(/)| 4 > |v|'I , |v/ , for 4>, 4/ E V 7 (4.12) 

for i = 1,2 and <7, <7 E Q- Finally, assume that 

<7 h-> (Bu + F)(t;q) is continuous from Q to L 2 ((0,7’), V 7 *) . (4.13) 

For arbitrary q N such that q N — ► q in Q, one then has the convergence 

z N (t\q N ) — > z (t;q) in V 7 norm 

, , „ (4.14) 

z?(t.;q N ) — > z t (t\q) in V norm 

for t > 0. Here 2 and z t are solutions to the linearized system (2.19) and z N and z^ 
solve the corresponding linear finite dimensional system in H N . 
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In our problems the admissible parameter, spare Q is taken to be a compart subset of 
the metric space Q — [L fXJ (0,a)] 3 x [III s ] 2 with elements satisfying the coefficient constraints 
specified in (4.7) as well as the physical constraints p > 0 , El > 0 and cpI > 0 on (0, a). In 
considering the remaining hypotheses, we note that the sesquilinear forms satisfy the ellipticity, 
continuity and parameter continuity conditions as long as damping is assumed in some region 
Cl C ft (see (2.14) and (2.14)). Moreover, the input term (E + Bu)(q) satisfies the condition 
(4.14). Finally, the convergence condition specified in (4.11) is satisfied as a consequence of 
the approximating properties of the cubic splines and Legendre polynomials in a (Jalerkin 
setting (see [6] for further details). 

Hence, for the' linearized problem with acoustic damping in D, a subsequence of solutions 
<y jN to the problems involving the minimization of the functionals (4.8) with v — 0. 1 , (4.9) or 
(4.10) subject to (4.2) will converge to a solution q of the original problem of minimizing the 
functionals (4.4), (4.5) or (4.6) subject to (2.19). The convergence in the case' involving the 
minimization of (4.8) with acceleration data does not follow directly from this theorem but 
can be obtained using results from the theory of analytic semigroups in a manner analogous 
to that used in [7]. 

With boundedness and Lipschitz continuity assumptions on the nonlinear coupling term 
(j(Z), similar results can be obtained for the nonlinear problem as summarized in the following 
remark. 


Lemma 2. (Nonlinear Damped System) Consider the system with the nonlinear input 
term Bu(t) + E(t) + G{Z(t)). If, in addition to assuming a continuity condition of the 
form (4.14) (with Bu + F replaced by Bu + F + (1(2)), we also assume that G{Z{t)) is 
continuous in t, uniformly Lipschitz continuous in 2, and displays at most, affine growth 
at oo, then convergence results analogous to those summarized in Theorem 2 can be 
obtained for the nonlinear system. 

We point out that these continuity assumptions on G were also made when discussing the well- 
posedness of the nonlinear system. Details concerning these conditions as well as arguments 
leading to the proof for the nonlinear case can be found in [5]. 

As indicated previously, the assumption of medium damping inside the cavity is often 
inappropriate in applications of interest. While we have not extended Theorem 2 to include 
the case in which acoustic damping is omitted, extensive numerical tests have indicated that 
parameter convergence and continuous dependence of the parameters on data is being obtained 
in the same manner exhibited by the system having both acoustic and structural damping. 
This is demonstrated by the results in the following examples for the nonlinear 2-D structural 
acoustic system in which cavity damping is omitted. 
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5 Numerical Examples 

To test the parameter estimation methodology, the general problem 
< i>tt — c 2 A(j) (x, y) £ Q , t > 0 , 

V0 • h = 0 (x, y) G F , /. > 0 , 


V<?i>(f,x,u>(f,x))- n = pfio t (t,x) 0 < x < .6 , / > 0 , 


pwtt + 


d 2 


dx 


,.2 


El 


d 2 w 


dx 


r.2 


+ CD I 


d^w 


dx 2 dt / 


4 


= — 4>t(t, X, w(t., x)) + f(t,x) + — - y] K.f Uj(t)Xpe,( x ) 

i = 1 

tn(<,°) = ^(/,°) = M^.6) = ^(<,.6) = 0 <>0 , 

</(0, x, y) = ^ ( (0,x, jy) = t«(0,.r) = re((0,x) = 0 


0 < .r < .6 , 
t > 0 , 


was considered. The cavity was assumed to have x and y dimensions a — .6 m and ( — 1 ill 
with a beam at one end having length .6 m, width .1 m and thickness .005 m (see Figure 3). 
The density and Young’s modulus for regions of the beam devoid of patches were taken to be 
p b = 2700 kg/m 3 and E = 7.1 x lO 10 N/m 2 which yields p = 1.35 kg/m and El = 73.96 N m 2 
for the linear mass density and stiffness parameter (see Table 1 for a compilation of the 
structural parameters for the system). The Kelvin- Voigt damping parameter was chosen to 
be c D I = .001 kg m 3 / .sec. Finally, the values pj = 1.21 kg/m 3 and c = 343 m/sec. were used 
for the atmospheric density and speed of sound. 

In the examples, we consider a system in which the bounding end beam has bonded 
to it a centered piezoceramic. patch covering 1/3 of its length as shown in Figure 3. The 
patch is assumed to have thickness T = .000508 m and width b = .1 ra (we point out 
that the chosen thickness value corresponds to 20 mil which is a commercially available 
thickness for piezoceramic patches). The Young’s modulus and density were taken to be 
E pe = 6.3 x lO 10 N/m 2 and p pc = 7650 kg/m 3 which are reasonable for a patch made from 
(1-1195 piezoceramic material. 

From (2.7) and (2.8), we see that the density and stiffness coefficient in the region of the 
combined beam and patch (Region 2) will be greater than that of the beam (Region 1) (see 
Figure 3). We also assume that the damping coefficient will be slightly larger in Region 2 
than Region 1 . 

For testing purposes, the structural parameter values in regions covered by the patches were 
chosen as specified in Table 1. As seen there, the constant K' B = E pe bd^(h + 27^ + T), which 
arises when modeling the actuation due to the patch, was taken to be K' B = .0067 Nm/V 
(this latter value was obtained by assuming a bonding layer of thickness T b e = .0001 m and 
taking c/ 31 = 1.9 x 10 -10 rn/V which is the value specified for (1-1195 piezoceramic material). 
The constant K' s which arises when using the patches as sensors was taken to have the value 
K' s = 170 V. 
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1 


Region I 
Region 2 

o 0.2 0.4 °- 6 

Fig. 3. Acoustic chamber ' with one centered 1/3 length patch . 

The following examples demonstrate the numerical estimation of the material parameters 
p, El ,ci)I ,K B and K s using various techniques for exciting the system and observing the 
response. In the first example, the natural frequencies for the fully coupled system were de- 
termined by simulating an impulse hammer impact to the center of the beam. The knowledge 
of these natural frequencies was then used when choosing the frequencies of the exciting forces 
in the other examples. In the second example, a periodic uniform forcing function (model- 
ing a uniform exterior acoustic pressure field) was applied to the beam for a short interval 
of time and then set to zero. This forcing function was chosen so that three system modes 
were initially excited and then allowed to begin decaying due to the damping the beam. I he 
acceleration of the center of the beam was used as data for estimating the material parameters 
p, El and col (since no voltage was applied in this example and the patch was not used for 
sensing, K B and were not estimated). In the third example, the system was excited by 
the application of a periodic voltage into the patch. Again, the system was excited for a short 
time interval and then allowed to freely decay in energy. Acceleration data at the center of 
the beam was used to estimate the four parameters p , El ,cpl and K B . The patch was used 
both for actuating and sensing in the fourth and fifth examples. In Example 4, a periodic 
voltage was applied for an initial time interval after which the system energy was allowed to 
decay. During the decay interval, two sets of data were calculated and a comparison was made 
between the results obtained when each was used for recovering the parameters p, E I ,cpl , K 
and K' s . The first data set consisted solely of the voltage produced by the patches during vi- 
bration (and hence contained strain measurements) while the second contained a combination 
of voltage measurements from the patch and acoustic pressure values from inside the cavity. 
Finally, a simulated voltage spike to the patches was used to excite the system in Example 5 
(with an effect similar to that observed when an impulse hammer is used to excite the system) 
with the patches again being used as a sensor throughout the remainder of the time interval. 
Thus in the last two examples, the “smart material” aspects of the structure were utilized in 
determining its physical parameters. 







True Values 

Initial Values 

Initial Values 




(Set 1) 

(Set 2) 

p 

Region 1 

1.35 

1.4 

1.4 

(kg/m) 

Region 2 

2.115 

2.0 

2.0 

El 

Region 1 

73.96 

74.0 

75.0 

( Nrn 2 ) 

Region 2 

125.4 

125.0 

127.0 

CpI 

Region 1 

.001 

.0008 

.0001 

(Nm 2 sec) 

Region 2 

.00125 

.0008 

.0001 

IC B (N 

m/V) 

.0067 

.007 

.01 

>c s (V) 

170.0 

172.0 

175.0 


Table 1 . (a) True values of the material parameters; (b) Set 1 of initial guesses for the material 
parameters; (c) Set 2 of initial guesses for the material parameters . 


Example 1: System Dynamics 

In order to determine the system dynamics with the parameter values in Table 1, the 
forcing function / was chosen to simulate an impulse at the center of the beam; that is, 

f{t, x) = 8(x - .3)8(t) . 

This models the force that would be delivered by a centered impulse hammer hit. After the 
initial impulse, the system was allowed to run unforced through time T — 8/60. 

The beam acceleration obtained with m x = m y = 12 and n = 16 basis functions at the 
point X — .3 is plotted in Figure 4 with a corresponding frequency plot in Figure 5. The first 
four system responses occur at 62.9, 179.7,342.1 and 397.9 hertz. 

For comparison, we note that the analytic natural frequencies of the first two symmetric 
modes of an isolated, homogeneous (no patches), undamped beam having the same dimensions 
as those in this system are 73.2 and 395.6 hertz while those of the uncoupled acoustic cavity 
are 171.5 and 343 hertz (the analytic undamped beam and cavity frequencies are given by 



respectively). The differences between the observed and analytic values are due to the presence 
of the patches on the beam, the internal damping in the beam, and the coupling between the 
beam and acoustic field dynamics. 
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The frequencies observed here can also be compared with the system values 65.9, 181.3, 
343.9 and 387.8 hertz which were obtained when no patches were bonded to the plate and 
linearized coupling conditions were assumed (see [8]). It can be seen that the increased stiffness 
due to the presence of the patches manifests itself at the higher frequencies (397.9 hertz versus 
387.8 hertz for the uniform beam) whereas the increased density is more of a factor at the 
lower frequencies (62.9 hertz versus 65.9 hertz for the uniform beam). This can be explained 
by comparing the bending shapes of the first and third beam modes and noting that the 
increased stiffness is more of an influence on the third beam mode while the increased density 
more directly influences the first mode. 


Beam Response to an Initial Impulse 



Fig. 4. The beam response to a centered impact. 


Beam Response to an Initial Impulse 



Fig. 5. The beam frequency response to a centered impact. 
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Example 2: Periodic Acoustic Excitation, Acceleration Data 


In this example, a uniform (in space) periodic force modeling an acoustic plane wave was 
used to excite the system and acceleration data was used to estimate the parameters p, El 
and c D l ( k' B and K' s were not estimated since we were not applying voltage to the patch or 
using the patch for sensing). Specifically, the forcing function was taken to be 


sin(1207rf) + sin(3607rf) -f sin(8007r<) , 0 <1/60 

0 , 1/60 <t. <8/60 


which initially excites the first, second and fourth system modes and then allows the oscillations 
to begin dying away due to the damping in the beam (see Figures 6 and 7). 

The parameter estimation was performed with data which was generated by calculating the 
acceleration of the central point of the beam at 498 uniformly distributed points throughout 
the time interval [0,8/60] (hence u = 2 in (4.8)). The acceleration was determined by using 
a second-order central difference on the displacements which were obtained by solving the 
nonlinear finite dimensional system (3.1) and evaluating w N at the points (.3, <*,), t ^ = 
k • — , k = 2, ■ - ■ , 499. Due to the relatively small number of frequencies being matched, all 

identification procedures were performed in the time domain which implies that C 2 in (4.3) 
was taken to be the identity. 

To test the algorithm and software, noisefree data was first generated using 120 acoustic 
and 15 beam basis functions (m x = m y = 10, n = 16). Using the initial parameter choices 
in column 5 of Table 1, the parameters were estimated using the same basis choice with the 
optimization being performed via a Levenberg-Marquardt routine. As demonstrated by the 
recovered values in the third column of Table 2, very accurate estimates of the parameters 
can be obtained when using the same number of basis functions when generating the data 
and estimating the physical parameters. 

To provide more realistic simulations, data was then generated using 168 acoustic and 
19 beam basis functions (m x = m y = 12, n = 20) with 120 acoustic and 15 beam basis 
functions again being used for estimating the parameters (the use of a larger number of basis 
functions when generating the data has the effect of adding numerical noise to the values being 
approximated in the optimization routine). The parameter estimates obtained with noisefree 
data and data to which 10% relative noise was added are reported in Table 2. While these 
results were obtained with the initial values in column 5 of Table 1, they are representative 
of those obtained with a variety of initial guesses (for this method of system excitation and 
observation, the optimization routine converged for a relatively large range of initial values). 
Time and frequency domain plots of the data containing 10% noise and the acceleration 
obtained with the estimated parameter values of Table 2, column 5, are given in Figure 6 
and 7. We note that in these figures it is essentially impossible to distinguish between the 
time data and the model response with estimated parameters with the graphical resolution 
used. These plots reinforce the observation that reasonable estimates of the parameters can be 
obtained with acoustic excitation and acceleration data with or without noise in the simulated 
data. 
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No Noise 
(Data Set 1 ) 

No Noise 
(Data Set 2) 

10% Noise 
(Data Set 2) 

p 

Region 1 

1 .3500 

1 .3199 

1.3469 

(kfi/m) 

Region 2 

2.1150 

2.1149 

2.1 181 

El 

Region 1 

73.960 

73.836 

72.130 

(Nth 2 ) 

Region 2 

125.100 

125.584 

130.716 

CD I 

Region 1 

.001000 

.001005 

.001023 

(Nin 2 scc) 

Region 2 

.001250 

.001232 

.001091 


Table 2. Estimated values of p, El and epl; (a) Data generated with m x = rn y = 1 0, n = 16 and 
no noise added to data; (b) Data generated with m x — m y — 12, n — 20 and no noise added to data; 
(e) Data generated with m x — m y — 12,7/ = 20 and 10% noise added to the data. 


System Response to an Initial Periodic Force 



Fl(i. 6. The system response to an initial periodic force. The solid line denotes the data containing 
10% nois( while the dashed line is the model response obtained with the recovered parameter values 
in column 5 of Table 2. 




System Response to an Initial Periodic Force 



Fig. 7. The system frequency response to an initial periodic force. The solid line (with x ’$) denotes 
the data containing 10 % noise while the dashed line (with o’s) is the model response obtained with 
the recovered parameter values in column 5 of Table 2 . 
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Example 3: Periodic Patch Excitation, Acceleration Data 

A second means of exciting the system is through the application of an out-of-phase voltage 
to the piezoceramic patches. The driving voltage* here was taken to be 

f si 11 ( I207t/) + sin(3607r/) + (.1 ) sin(800^/) , Q<t<l/6Q 

u(t) = < 

[ 0 , 1/60 < t < 8/60 

which initially excites the first four system modes and then allows the oscillations to begin 
dying away due to the damping in the beam (see Figures 8 and 9). Because the patch excites 
the* higher frequencies more efficiently than did the simulated sound field of the last example, 
the magnitude of the 400 hertz contribution was reduced to better balance the system response. 

Acceleration data obtained in the same manner described in the last example was used to 
estimate the four parameters p, EI^cpI and K' B and the estimates obtained with no noise and 
10% noise added to the data are summarized in Table 3 (see Table 1 for the true parameter 
values). Figures 8 and 9 contain the time and frequency domain results for the case in 
which the data, contains 10% noise. In both cases, the choices m x = m y — 12, n — 20 
and m x = rn y = 10, n = 16 were used when generating the data and estimating parameters, 
respectively. Also, while similar results were obtained with a relatively large variety of starting 
parameter values, the parameter values reported here were obtained with the initial choices 
in column 3 of Table 1 . 



No Noise 

10% Noise 

p 

Region 1 

1.3502 

1.3425 

(ky/m) 

Region 2 

2.1211 

2.0973 

El 

Region 1 

73.139 

73.449 

(Nm 2 ) 

Region 2 

128.179 

124.741 

C D i 

Region 1 

.000997 

.000849 

( N in 2 sec) 

Region 2 

.001256 

.001787 

K B (N 

rn/V) 

.006771 

.006654 


Table 3. Estimated values of p, El, C[jl and (a) No noise added to the data; ( b ) \0% noise 
added to the data . 


In comparing the results obtained with noisefree data with those obtained under the same 
conditions (same number of basis functions and noisefree data) in the last example, we note 
an approximately 2%! change in the estimated value of the stiffness parameter in the region 
covered by the patches. This difference appears to be due to the fact that the “true” data 
is calculated with a larger number of basis functions and hence with more accuracy than the 
solutions being obtained during the parameter estimation. The resulting “numerical noise" 
manifests itself more strongly in this case since the patches are more effective than the acoustic 
field at exciting high frequency oscillations which require greater accuracy to resolve (this 
tendency is also noted in tin* remaining examples where patch activation is used to excite the 
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system). We point out that when the same number of basis functions are used for generating 
data and estimating parameters, the results are essentially identical to those in column 2 of 
Table 2 with highly accurate estimates of the physical parameters. Finally we note that the 
“numerical noise” due to the differing discretizations combines with random noise added to the 
data to simulate the noise which is present when real data is used to estimate the parameters 
in a physical experiment. 


System Response to an Initial Periodic Voltage 



Fig. 8. The system response to an initial periodic voltage. The solid line denotes the data containing 
10% noise while the dashed line is the model response obtained with the recovered parameter values 
in Table 3. 


System Response to an Initial Periodic Voltage 



Fig. 9. The system frequency response to an initial periodic voltage. The solid line (with x s) 
denotes the data containing 10% noise, while the dashed line (with o’s ) is the model response obtained 
with the recovered parameter values in Table 3. 
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Example 4: Periodic Patch Excitation, Patch and Pressure Data 


In the previous two examples, acceleration data modeling that which would be obtained 
from a centered accelerometer was used in the estimation of the physical parameters. In this 
as well as the next example, both the excitation of the system and the sensing of the beam 
dynamics are performed with the patches, thus utilizing their “smart material 11 capabilities. 
The beam data in this case consists of the voltage that is produced when the patches strain 
during vibration (hence the criterion functional (4.9)) as well as data consisting of voltage 
measurements in combination with pressure measurements from within the cavity (with the 
criterion functional (4.10)). 

Here a periodic driving voltage was applied to the patches for the first 1/5 of the time 
interval [0, 10/60] of interest after which the system was allowed to begin decaying in energy. 
Specifically, the input voltage was taken to be 

f si n ( 1 207 t/ ) + sin(3607r/) -f sin(8007r/) , 0 < / < 2/60 

u(t) = < 

{ 0 , 2/60 < t < 10/60 

with the frequencies again chosen so as to strongly excite the first, second, fourth and sixth 
system modes (see Figure 1 1). During the final 4/4 of the time interval, the voltage produced 
by the patches as well as the acoustic pressure at the cavity point (.6,.l) were calculated at 
750 uniformly distributed times throughout the interval [5/120,10/60]. These values were 
then used as our simulated data. 

The estimated values of the parameters parameters p, E I , cdI , K B and K s are recorded in 
Table 4. In the first four simulations, 168 cavity and 19 beam basis functions were used to 
calculate the data while 120 cavity and 15 beam basis functions were used in the estimation 
of the physical parameters. The values in columns 5 and 4 were obtained with data, consisting 
of the 750 voltage values generated by the piezoceramic patches while data consisting of both 
voltage and acoustic pressure values was used to obtain the results in columns 5 and 6. Finally, 
the results from a simulation in which 120 cavity and 15 beam basis functions were used for 
both the generation of data and estimation of parameters are reported in column 7 of Table 4. 
The data for this latter simulation consisted of the previously described combination of voltage 1 
and pressure values. 

It is first noted that when the same number of basis functions are used for tlu 1 generation 
of data and estimation of parameters, highly accurate results can be obtained for a, variety of 
initial starting values when using a combination of voltage and pressure values as data (the 
results in column 7 of Table 4 were obtained with the initial values in column 5, Table 1). 
This is consistent with the results reported in Example 2 and described in the Example 3. As 
discussed in the latter example, the use of a larger number of basis functions when determining 
the data adds a form of “numerical noise 1 ’ since the data is calculated with greater accuracy 
than are the solutions obtained during parameter estimation. Hence, although no random 
noise was added to the data, the estimated values in columns 3,5 and 6 of Table 4 differ 
slightly from the “true” values listed in Table 1 but are consistent with those obtained in the 
last example (see column 3, Table 3). 

The effects that the accuracy of the initial guesses had on the estimated parameters can 
be seen by comparing the results in columns 3 and 4 (voltage data) as well as in columns 5 
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and 6 (voltage and pressure data). The parameter values in columns 3 and 5 were obtained 
with the highly accurate initial choices of column 4, Table 1 while those in columns 4 and 
6 were calculated with the less accurate values in column 5, Table 1 which could reasonably 
have been obtained by first visually fitting the data in the frequency domain. As observed, the 
change in starting values led to fairly large changes in the recovered parameter values when 
purely voltage data was used, whereas almost no change was noted when the data consisted 
of voltage and pressure values. 

The variation in the recovered parameter values obtained with voltage data can partially 
be attributed to the flatness of the criterion functional with strain observations as compared 
to that seen with acceleration data. The inclusion of pressure values in the data adds rich- 
ness as a result of the added information about the acoustic state as well as the fact that 
higher frequencies are more easily observed in the pressure data than in the stain values (see 
Figure 10). Due to the added information in the data containing strain and pressure values, 
a wider range of initial guesses could be used since the optimization routine was less like to 
become stuck in local minima. 



Sim 1 

Sim2 

Si m3 

Sim4 

Si in 5 

p 

Region 

1 

1.3501 

1.3745 

1.3502 

1.3502 

1 .3500 

(kg/jn) 

Region 

2 

2. 1207 

2.1476 

2.1230 

2.1231 

2.1150 

El 

Region 

1 

73.178 

76.319 

72.910 

72.908 

73.960 

( Nm 2 ) 

Region 

2 

128.022 

123.546 

129.041 

129.051 

125.400 

CD I 

Region 

i 

.000992 

.000283 

.000989 

.000988 

.000999 

(Nrn^sec) 

Region 

2 

.001274 

.000431 

.001295 

.001297 

.001250 

K' B ( N 

m/V) 


.007274 

.006709 

.006767 

.006766 

.006700 

X s (V) 

158.749 

85.832 

172.018 

171.977 

169.999 


Table 4. Estimated values of p , E I ,cpl , IC B and IC^; (a) Sim l : voltage data and the initial guesses 
of Table 1, column 5; (b) Sim2: voltage, data arid the initial guesses of Table 1, column 6; (<:) SimS: 
voltage and pressure data with the initial guesses of Table 1 , column 5; (d) Sim 4: voltage and pressure 
data with the initial guesses of Table l } column 6; (e) SimS: voltage and pressure data generate with 
same number of basis functions used in parameter estimation , initial guesses of Table /, column 6 . 
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System Response to an Initial Voltage - Beam Strain System Response to an Initial Voltage - Acoustic Pressure 





System Response to an Initial Voltage - Acoustic Pressure 



Fig. 10. The system response to an initial periodic voltage. The solid line (with x's) denotes the 
data containing 10% noise while the dashed line (with o’s) is the model response obtained with the 
recovered parameter values in Table J h column 6 . 


Example 5: Patch Impulse Excitation, Patch Data 

The input voltage in this example was taken to be a narrow triangle which simulates a 
voltage impulse to the patches. As demonstrated by the strain and pressure plots in Figure 12, 
this causes the excitation of multiple system modes which through time decay in energy due 
to the damping in the beam. 

The data was generated using 168 cavity and 19 beam basis functions, and a comparison 
was again made between the results obtained with data consisting solely of voltage values 
(using the criterion functional in (4.9)) and data made up of both voltage and pressure mea- 
surements (with the criterion functional in (4.10)). In both cases, the data was calculated at 
900 uniformly spaced points throughout the time interval [1/60, 10/60] 

The parameter values obtained with 120 cavity and 15 beam basis functions and strain 
(voltage values) data are given in Table 5. The initial guesses of column 5, Table 1 were 
used to obtain both the results containing 10% random noise and those to which no noise was 
added (the noisefree results can be seen to be quite close to those obtained in the last two 




examples). Finally, plots containing the time and frequency domain results for the case in 
which 10% random noise was added to the data are given in Figure 12. 

We point out that with this means of system excitation, a much wider range of initial 
guesses led to convergence to the reported values than was the case in the last example where 
a periodic voltage was used to excite the system. This appears to be due to the fact that a 
larger number of frequencies are excited and observed in the strain data which in turn leads 
to “richer” data. On the other hand, the large number of frequencies observed in the pressure 
field ( see Figure 12) made it difficult to use pressure values in combination with the voltage 
as data since an extremely good initial guess was required in order to obtain convergence 
of the optimization routine (in this case, one must work nearly exclusively in the frequency 
domain at first since any deviation in the density and stiffness parameters leads to frequency 
changes (see (5.1)) that make estimation in the time domain very difficult). We also noted 
that similar problems were encountered when acceleration data was used with this number of 
frequencies excited due to the fact that the acceleration provides a more sensitive measure of 
beam movement than do the patches which measure accumulated strain. Hence, with a large 
number of frequencies excited, the role of the patches as the sole sensors in the structural 
acoustic system appears to improve, and in some cases this method of data observation may 
be preferable since the matching of a smaller number of observed frequencies may lead to a 
more tractable optimization problem. 




No Noise 10% Noise 

p 

Region 1 

1.3516 

1.3534 

(kg/m) 

Region 2 

2.1236 

2.1218 

El 

Region 1 

73.261 

73.969 

( Nm 2 ) 

Region 2 

128.156 

126.152 

Cl) I 

Region 1 

.000993 

.000959 

( Ntn 2 sec ) 

Region 2 

.001273 

.001334 

k: b {n 

mi V) 

.006677 

.006188 

k: s (V) 

173.200 

183.266 


Table 5. Estimated values of p, E I , cpI , K B and /C s using strain data; (a) No noise added to the 
data; (b) 10% noise added to the data. 
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Linear Amplitude 






6 Conclusions and Implementation Issues 

In this paper, a parameter estimation methodology for structural acoustic systems in which 
piezoceramic. patches on the structure act as sensors and actuators has been presented. While 
illustrated throughout this presentation in the context of a 2-D acoustic cavity with a thin 
beam at one end, the methods developed here are equally applicable in many 3-D distributed 
parameter models (such as those in [9, 10]) of actual experimental devices currently being used 
to test various modeling, parameter estimation and control schemes involving piezoceramic 
actuators and sensors. Recently, these estimation methods have been successfully used as an 
integral part of feedback control techniques with experimental data - these findings will be 
reported in a future report. 

The emphasis in developing this parameter estimation methodology was on providing a 
scheme which was amenable to approximation and implementation under a variety of damping 
assumptions. In doing so, conditions leading to well-posedness and parameter convergence 
results for both the linearized and original nonlinear problem were formulated, and suitable 
numerical techniques for approximating system dynamics and implementing the parameter 
estimation schemes were developed. 

To illustrate the method, several numerical examples illustrating a variety of techniques for 
exciting the system and generating data were presented. Throughout the examples, emphasis 
was placed on simulating exterior forces and generating data in a manner consistent with that 
used in actual structural acoustics experiments and applications. This was done so as to gain 
insight regarding the effectiveness of the parameter estimation method in various settings. The 
first step in the examples was the determination of the natural frequencies for the symmetric 
(in x) modes of the coupled system. This knowledge was then used when determining driving 
frequencies so as to evoke particular system responses in the later examples. Moreover, by 
comparing the results for the system under investigation with comparable results obtained for 
the structural acoustic system in which no patches were bonded to the beam as well as analytic 
values of the natural frequencies of the uncoupled beam and acoustic wave, the qualitative 
effects of coupling, damping, and material change's due to the presence of the patches were 
determined. 

A simulated acoustic source was used to excite the beam, and hence the system, in the 
second example with data consisting of acceleration values calculated at the center of the 
beam. While the numerical simulations demonstrated the success of the method for both 
noisefree data and data to which random noise had been added, this means of exciting the 
system in a manner that can be accurately simulated will be difficult to implement. This 
is due to the fact that although speakers can be used to create multiple frequency acoustic 
forces of the type used in the example, it will be difficult, if not impossible, to cleanly cut the 
acoustic excitation at a given time since the speaker and room will continue to reverberate 
and echo even after the pow r er is cut. This echo or reverberation will be difficult to simulate in 
the system model thus making parameter estimation through this means of system excitation 
difficult to use in practice. 

A more readily simulated means of exciting the system is through the application of a 
prescribed voltage to the patches and this was the input in the remaining examples (by altering 
the frequencies and magnitude, this voltage can be tailored so as to evoke a desired system 
response). In the third and fourth examples, a numerically simulated multiple frequency, 
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periodic voltage was used to initially excite the system at which point the voltage was cut and 
the system energy was allowed to begin decaying due to the damping in the beam. In addition 
to being more accurately simulated numerically, this also proved to be a more effective means 
of exciting high frequency responses in the system than was the acoustic excitation. A short 
duration triangular input was applied to the patch in the final example to simulate an initial 
voltage spike. 

In the third example, acceleration values on the beam wore used as data, whereas voltage 
values (measuring accumulated strain in the beam) and voltage values in combination with 
acoustic pressure values were used as data in the fourth and fifth examples. In these examples., 
it was found that when only a few system frequencies were excited (three or four), the accel- 
eration data was sufficiently sensitive so as to permit effective parameter estimation whereas 
highly accurate initial guesses wen* needed in order to use beam data containing only strain 
information (which qualitatively has the properties of displacement data). By augmenting the 
strain data with pressure measurements, however, sufficient information was added so as to 
again lead to successful parameter identification with a range of initial values. On the other 
hand, the excitation of large 1 number of system frequencies through the* simulated voltage' 
spike to the* patch led to strain data which was sufficiently “rich” so as to permit success 
of the method while* the introduction of pressure data led to a failure in the optimization 
routine due to the very large number of frequencies in the pressure* measurements. liene e the 
success of the* method with strain measurements comprising part of the data depended partly 
on the number of excited frequencies; for a small number, the strain measurements had to be 
supplemented with acceleration or pressure values in order to ensure parameter convergence 
for data containing any noise. 
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